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Abstract 

The coherent interband dynamics of optically excited two-dimensional lat- 
eral surface superlattices driven by an in-plane static electric field has been 
investigated. The linear absorption, the spectrally-resolved pump-probe four- 
wave mixing signals and spatial coherent wavepacket evolution in the time- 
domain are obtained. When the rational condition E x /E y = p/q, with p, q 
prime to each other, is fullfilled, it is found that p peaks appear within the 
frequency interval Ub x = eE x a/U in both linear absorption and degenerate 
four-wave mixing signals. The coherent time evolution of the electron-hole 
pair wavepacket is that of a breathing mode, with the period 2ttp/lobx, These 
findings are consistent with the recent spectral results (Phys. Rev. Lett. 
86, 3116), hence providing a method for probing the coherent dynamics of 

quantum particles in 2D lattices. 

PACS numbers: 78.47.+p, 73.50.Fq, 73.20.Dx 
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I. INTRODUCTION 



Recently, there have been much investigation on the electronic and transport properties 
of semiconductor super lattices driven by the external static and/or time-dependent electric 
fields Jl [14|]. It is foreseen that these explorations will uncover rich physics and even find 
potential applications in electro-optic devices. In this context, a lot of interesting theoretical 
and experimental findings have been reported. Among these, we mention Bloch oscillations 
M, Fano resonances 0, absolute negative conductivity ||, inverse Bloch oscillators [0, 
dynamic localization MM, dynamic Franz-Keldysh effect M, self- induced Shapiro effect JTU 



dynamic fractional Wannier-Stark ladders [fLlJ| . Particularly, with the development of atom 
optics, physicists began to link solid state physics and optics together, both theoretically 
and experimentally Hl5| . The recent clean demonstration of Wannier-Stark ladders 0] and 
Bloch oscillations |l7j in the atom optics are typical examples. 

Several months ago, Gliick et al calculated numerically the ground Wannier-Bloch band 



of a quantum particle moving in a lattice under an in-plane electric field |p.8| . The lattice 
may be an array of quantum dots for electrons, or an interference pattern of light waves for 



cold atoms ||15|| . They found that the energy spectrum resembles the Hofstader butterfly 
||T9| in a magnetic system, and the 1-D superlattice driven by dc-ac electric fields [[LI]]. They 
showed that for a rational static force F of components in the ratio F x /F y = p/q, with p, q 
coprime integers, there are s = p 2 + q 2 identical subbands separated by the energy interval 
5 = \F\a/ y^s, where a is the lattice constant. Since the lattice period in the direction of the 
field is \fsa, this result predicts a splitting of the 1-D Stark period into s equally spaced 
subbands. The question arises as to whether such a spectrum is measurable. 

In this paper we report our finding that this striking splitting in the spectrum can be 
identified in optically excited two-dimensional lateral surface superlattices driven by an in- 
plane static electric field E. Such 2-D superlattices can be prepared by embedding a square 
array of GaAs cylinders in a very thin Ga^Ali-^As epilayer p0| . By numerically solving 
the generalized semiconductor Bloch equations, this spectrum is identified in both the linear 
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absorption signal and spectrally resolved pump-probe degenerate four-wave mixing signals. 
Also, we link the breathing-mode period of spatial coherent wavepackets oscillations in the 
time domain in the system to the splitting of the spectrum. 



II. MODEL AND NUMERICAL RESULTS 



We consider an electron system moving in a 2-D lattice and an external in-plane electric 
field E with its component in the ratio E x /E y = p/q, where p and q are coprime numbers. 
The energy dispersion is treated in a two-band model, the valence and conduction bands. 
We discuss the dynamics in the context of the following semiconductor Bloch equations 
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where Pk(t) is the interband polarization and n e (M k(t) the electron (hole) population density 



in the conduction (valence) band. The quantities e iik = e^k — Dq^jk-qi^ 



e, h) 



are the renormalized electron and hole energies due to the Coulomb interaction. Also, 
^R,k = (d cv £ + X)k M k -q|-^q)/^ are ^ ne renormalized Rabi frequencies, with d cv the dipole 
moment and S(t), the Gaussian laser pulses. The relaxation time approximation has been 
assumed, with and Tj- being the longitudinal and transverse relaxation rates, respectively. 

The spectrally-resolved absorption, pump-probe nonlinear four-wave mixing signals and 
the real-time coherent evolution of the wavepackets in real-space have been calculated nu- 
merically by solving the coupled Eqs.(l) and (2). The energy bands of 2-D lateral sur- 



face superlattices of square geometry can be described by the tight-binding form |20| . [21| . [23 



Cj(k) = -^(cos(k x a) + cos(k y a)) (i = e, h), with A C (A„) being the combined miniband width 
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of the conduction (valence) band in the x and y directions. For simplicity, an on-site Coulomb 
interaction has been adopted, which describes the first-order Born scattering among the car- 
riers. Although the excitonic interaction is simplified to be on-site in this work, it is expected 
that the inclusion of more realistic models will not change our findings qualitatively. This 
kind of approximation has been used with success by many groups [^-[4] |Hi|fMH . 

Linear absorption: In the linear response regime, one can completely neglect Eq.(2), and 
set the electron (hole) density n c („)(k, t) to vanish in Eq. (1). The absorption a(uj) is pro- 
portional to the imaginary part of the first-order susceptibility X«M=Im[p( 1 )H/E(a;)], 
where PW(w) is the Fourier transform of the quantity J2k P (k, t). The integro-differential 
equations can be solved numerically in the accelerated basis k — | J * E(r)c/r. Parameter 
values used in the calculations are A c = A„ = 20 meV, V = 10 meV, and T L = T T = 0.5 
THz. The full width at half maximum for the strength of the Gaussian laser pulse is as- 
sumed to be 59 femtoseconds. In the simulation the Bloch frequency ujb x = eE x a/h of the 
static electric field in the x direction remains unchanged and taken to be 37r THz, while 
the Bloch frequency u>b v = eE y /h in the y direction takes three different values: 1/3 u Bx , 
2/3 ojbxi and 4/3 ujbx- Distinct peaks are clearly identified in the three panels, where the 
ratio E x /E y is indicated in the right-up corner. Roughly three peaks appear within every 
interval uj Bx in the three panels, although the Bloch frequencies in the y direction differ by 
a factor of four in the top and bottom panels. This property is the direct reflection of the 
striking spectrum discussed by Gliick et al fll8| . Since u Bx = eE x a/h = (p|eE|/ v /s)(a//i), 



and the characteristic energy splitting 5 = \eE\a/h^, then 5/ujbx = 1/p, where for Fig.l, 
p = 3 in all three panels. The interesting and special case of E x : E y = 3 : 3 = 1 : 1 is shown 
in the Fig. 2, where peaks appear only once in the fixed interval ujbx- It is clear from the 
comparison of Fig.l and Fig. 2 that the number of peaks appearing in the fixed frequency 
interval ujb x is solely determined by the integer p, instead of the absolute value of E x and 
E y . It should be noted that the absorption peaks obtained in the present model are much 
more stable that those in dc-ac field case [J25|J . The latter are easily blurred by the dephasing 
and Coulomb interaction- mediated Fano interferences ||. 



Four-wave mixing signals: Four-wave mixing experiments are usually applied to observe 
the coherent signals from optically excited semiconductor samples by taking advantage of the 
fact that this technique can substantially avoid the inhomogeneous broadening || . In order 
to mimic the experimental realization of probing the spectrum, the two laser pulses: the 
strong pump pulse propagating in the k x direction, and the weak probe pulse propagating 
in the direction k2, are arranged to be delayed a time r p relative to each other, 

£{t) = £ x (t) exp[i(ki ■ r - tot)] 

+ £ 2 {t - r p ) exp[z(k 2 • r - oo(t - r p ))] . (3) 



One can use the method of Lindberg et al [26], or that by Banyai et al [27|, to extract the 



nonlinear pump-probe four- wave mixing signals propagating in the direction 2k 2 — ki. The 
spectrally-resolved pump-probe degenerate four-wave mixing squared signals are shown in 
the four panels of Fig. 3 for different values of r p . A Gaussian profile for the pump and probe 
laser fields £\{t) and 82(f) is assumed, and the other parameters used are the same as those 
in Fig.l. In the figure, we only show the case of uj By /u Bx = 2/3. The other sample cases 
u>By/^Bx = 1/3,4/3 have similar shape, and are omitted for saving space. Inspection of the 
four panels clearly show that the peaks also appear three times in the frequency interval 
ujbx, although a different time delay t p is employed in different panels. This demonstrates 
that there exists the energy interval/ladders in the system whose energy spacing is identical 
and equal to u>b x /3. This is a direct manifestation of the split stable resonance spectrum 
found in 2-D lattices under an external in plane uniform field. 

Coherent wavepacket oscillations: As is well known, the time-domain wavepacket Bloch 
oscillations inferred from the semiclassical approach are a direct proof of the frequency- 
domain resonance states of Wannier-Stark ladders ]!]. While in Fig. 1-3 the signals are all 
from the frequency domains, what about the evolution of the coherent wavepacket in the 
time and spatial domain? In the following, we will show the coherent spatial wavepacket 
oscillatory behavior in the realm-time domain [ |30| . 

According to the definition of the electron-hole wavepackets |29 : P(r, t) = | J2u exp(— ik- 



r)i = i c (t) | 2 , where r is the relative position vector of the electron-hole pair. The coherent 
spatial wavepackets have been shown sequentially in Fig.4, with the corresponding time 
displayed in the upper-right corners. In the simulation we used the following parameters: 
^Bx = 67T THz, E y /E x = 2/3, while the other parameters used are the same as those in 
Fig.l. The contour of the wavepacket in the earliest time (t = 0.05 ps) looks like a small 
circular dot, which can be regarded as the profile of the initial excitation of the electron- 
hole pair. With the development of time (t = 0.20 ps), due to the driving electric field, the 
wavepacket is pulled apart and distributed along both the x and y directions in a anisotropic 
manner. When time reaches 0.5 ps, the wavepacket distribution along the y direction is 
reduced. This is expected because the Bloch oscillations period along the y direction is 
Tb v = 2'k/ujbx — 0.5 ps. For the same reason, the wavepacket distribution experiences 
narrowing along the direction x at time t = 0.67 ps(~ 2Tb x )- At time t = 0.8 ps, the 
wavepacket again inflates in both x and y directions. The most striking phenomena perhaps 
lies in the last panel, where the wavepacket has shrunk into a small constellate at time t = 1 
ps. Since T = 2irh/5 = 2ttp/ub x , for the frequency chosen one obtains T = 1 ps, which 
equals three times T Bx , and twice T By - This indicates that the energy splitting 5 found in 
Ref. [jTSJI is the direct consequence of the wavepacket Bloch oscillations with period T of 



the electron-hole pair, although no hole motions are involved in Ref. |[L8| . The existence of 
this period is thus a direct manifestation of commensurate synchronous Bloch oscillations in 
both x and y directions. It should be mentioned in passing that in the absence of excitonic 
interaction, the wavepacket will undergo the breathing mode in a precise time period T. 
The presence of the interaction changes only slightly this period. 

III. CONCLUDING REMARKS 

The linear absorption and nonlinear pump-probe spectrally resolved four-wave mixing 
signals have been calculated in optically excited 2-D lateral surface superlattices driven by 
an in-plane static electric field. The ratio of E x to E y is chosen to be the rational p/q, 
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making the peaks appear p times in the frequency interval u Bx, both in linear absorption 
and the four-wave mixing signal. This finding provides a way to test experimentally the 
spectrum arising in the dynamics of particles moving in a 2D lattice and a constant force 
||18|| . Also the coherent wavepacket oscillates with a period corresponding exactly to that 
associated with the energy interval 5, showing that the energy interval is a consequence of 
the commensurate synchronous Bloch oscillations in both the x and y directions. 



Theoretically, in the absence of all kinds of dephasing ||18|| , the spectrum has the property 
of extreme sensitivity to fluctuations in the direction of the field. For instance, if E y /E x = 
q/p, then an arbitrarily small rotation yields a ratio q'/p', with (q',p') — > oo. For instance, 
if q/p = 1/3, a tiny change in the field produces the fraction 1001/3000, while an even 
smaller change gives rise to the ratio 1000001/3000001, thus modifying wildly the splitting 
size as determined by the denominator fraction. The presence of various kinds of dephasing 
mechanisms blur this fine structure, making the identification of the peak in the fixed interval 
ujbx impossible, when p becomes large []3T|. Experimentally, the phenomena predicted here 
can be observed in a photoexcited 2D quantum dot array driven by an in-plane electric field, 
with the electric field strength tuned to be a rational. It is expected that clean observation 
of these ladders can be also realized in ultracold atoms dwelling in the two perpendicular 
optical waves, accelerated accordingly. 
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FIGURES 

FIG. 1. Linear absorption a(u>) plotted as a function of (u)—u) g )/ujBx, where uj g is the frequency 
associated with the conduction-valence band gap. 

FIG. 2. The same as Fig.2, except that us y ■ ujbx = 3:3 = 1:1. Only one peak appears in 
every fixed frequency interval ujb x . 

FIG. 3. Spectrally-resolved degenerate squared four- wave mixing signal \P^(lo)\ 2 plotted as 
a function of (uj — uj g )/ijjB,x- Three equidistant peaks appear within the fixed frequency interval 

UBx- 

FIG. 4. Coherent time-evolution of a wavepacket, showing the breathing behavior in the 
real-time domain. The oscillatory period is 1 ps, revealing the energy spacing in the fre- 
quency/energy domain. 
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